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NON-MAXWELLIAN  FREE-ENERGY  GENERATION  IN  THE  MAGNETOTAIL 
DUE  TO  CHAOTIC  PARTICLE  MOTION 


I.  Introduction 

The  notion  of  charged-particles  in  the  earth's  magr.etotail  has  Peer, 
under  extensive  investigation  for  the  past  two  decades  in  order  to  fully 
understand  the  properties  of  the  magnetotail  plasmas  from  both  the 
experimental  and  theoretical  points  of  view.  A  number  of  specialised 
aspects  of  the  single- particle  motion  have  Peer,  analyzed  using  various 
methods  including  approximate  analytical  methods  (Speiser,  1965;  Alexeev 
and  Kropotkin,  1970;  Sonr.erup,  1971;  Stern  and  Palmadesso,  1975;  Cowley, 
1973;  Peilat  and  Schmidt,  1979)  and  numerical  methods  (Speiser,  1967; 
Cowley,  1971;  Zastwood,  1972;  Swift,  1977;  Wagner  et  al . ,  1979;  Gray  and 
Lee,  1982;  Spei3er  and  Lyons,  1984). 


The  magnetotail  field  may  be  modeled,  in  its  simplest  form,  by  a 

A 

neutral  sheet  magnetic  profile  3Q(z)x  with  3Q(z  *  0)  -  0  and  a  superimposed 
component  3^z  normal  to  the  plane  of  the  neutral  sheet  (the  so-called 
''quasi-neutral  sheet"  geometry).  We  use  the  standard  magnetotail 
coordinate  system  with  :<  in  the  earthward  direction.  In  the  above  works, 
the  fundamental  question  regarding  the  ir.tegrability  of  the  particle  motion 
was  not  addressed.  However,  it  has  recently  been  shown  (Chen  and 
Palmadesso,  1984a)  that  the  magr.etotail-like  system  is  intrinsically 
nonir.tegrable  due  to  the  presence  of  the  normal  component  3.,. 
Subsequently,  Chen  and  Palmadesso  (1985),  henceforth  referred  to  as  C-P, 

have  identified  various  phase  space  structures  in  detail.  More 

0 

specifically,  U3ing  the  Poincare  surface  of  section  methoc,  they  have  shown 
that  particle  motions  in  the  magr.etotail-like  configuration  can  be 
classified  into  three  distinct  types  of  orbits  occupying  disjoint  regions 
of  the  phase  3pace;  the  bounded  ir.tegrable  orbits,  'unbounded  stochastic 
jroits  and  unbounded  transient  orbits. 

In  studying  individual  orbit3,  a  number  of  researchers  have  noted 
" randomness"  in  certain  orbit3.  For  example,  Swift  (1977)  noted  that  seme 
croits  appear  to  randomise  Lr.  the  equatorial  plane  after  a  few  crossings. 
Wagner  et  al.,  (’979)  found  that  certain  orbits  exhibit  sensitive 
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dependence  cn  initial  conditions.  Gray  and  Lee  (1982)  noted  that  the 
magnetic  moments  for  some  particles  exhibit  randomness  across  the 
equatorial  plane.  We  now  understand  that  this  randomness  is  a 

manifestation  of  the  intrinsic  r.onintegrable  nature  of  the  system.  The 
work  of  Chen  ar.d  Palmadesso  ( 1 98-^a;  *985)  provides  a  systematic  and  unified 
understanding  of  the  nature  and  properties  of  the  particle  motion 

throughout  the  entire  phase  space.  More  recently,  Martin  (1985)  showed 
that  particle  motion  near  an  X-poir.t  is  chaotic.  Also,  it  has  long  been 
known  (Dr ago  and  Finn,  1976)  that  motion  in  the  dipole  field  is  not 

integrable. 

The  integrability  of  the  particle  motion  is  not  purely  of  mathematical 
interest.  The  existence  of  distinct  classes  of  orbits  has  profound  effects 
on  the  plasma  particle  distribution  and  on  the  response  of  the  magr.etotail 
plasmas  to  external  influences.  In  thi3  regard,  the  concept  of 
"differential  memory"  was  introduced,  by  C-P  to  describe  the  property  that 
particles  in  the  disjoint  regions  of  the  phase  space  respond  to  external 
influences  or.  different  time  scales.  This  implies  that  a  plasma 
distribution  function  has  a  natural  tendency  to  develop  non-Maxwell i an 
features  in  response  to  changes  in  physical  conditions.  In  turn,  the 
distribution  and  character  of  the  charged-particle  orbits  car.  have  a 

significant  impact  on  other  dynamical  properties  of  the  sagnetotail.  One 
important  example  is  the  collisionless  tearing-mode  instability  (Furth, 
1962;  Pfirsch,  ’96 2;  Laval  et  al.,  1966)  which  has  long  been  thought  to 
play  an  important  role  ir.  magnetic  field  reconnection  (Coppi  et  al,  1966; 
Schindler,  1966).  It  has  recently  been  shown  that  the  growth  rates  of  the 
collisionless  tearing  instability  can  be  enhanced  by  up  to  a  few  orders  of 
magnitude  due  to  the  presence  of  temperature  anisotropy  or  other  non- 
Maxweiliar.  features  in  the  particle  distribution  (Chen  ar.d  Palmadesso, 
193^0;  Cher,  et  al . ,  198^;  Cher,  ar.d  Lee,  1985).  Thus,  the  particle  dynamics 
plays  e  fundamental  role  ir.  influencing  the  magr.etotail  dynamics.  A 
fundamental  question,  then,  13  whether  the  magr.etotail  car.  sustain  free- 
er.ergy  carrying  nor.— Maxwellian  features  that  car.  support  large-scale 
instabilities  of  potential  relevance  to  magr.etotail  dynamics,  such  as  the 
ncn-Maxwelliar.  collisionless  tearing  mode. 

In  this  paper,  we  will  examine  the  way  the  magnet. c  field  topology  of 
a  quasi-neutral  sheet  car.  generate  non-Maxwell:  an  features  in  plasma 
distribution  functions  and  the  attendant  free-er.ergy  that  car.  drive  certain 
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plasma  instabilities.  In  Sec.  II,  we  review  the  particle  dynamics  using 
the  Poincare  surface  of  section  method.  In  Secs.  Ill  and  IV,  physical 
implications  for  magnetcspheric  dynamics  will  be  discussed.  We  will 
primarily  refer  to  the  ion  motion  for  illustration.  The  modifications 
needed  for  electrons  are  trivial. 

II.  Single-Particle  Dynamics 

In  order  to  model  the  motion  of  charged  particles  in  the  magnetotail, 

A  A 

we  consider  a  magnetic  field  given  by  3  *  3Qf(z)x  *  3,,z,  where  30  is  the 
asymptotic  field  in  the  x  direction,  3n  is  the  uniform  normal  field  and 
30f(z)  is  a  neutral  sheet  profile  sucn  that  f(-z)  -  -f(z).  We  will 
primarily  use  the  Harris  configuration 

f (z )  -  tarJi(z/<S),  (*) 

where  <5  is  the  characteristic  scale  length  of  the  magnetic  field.  The 
treatment  i3  also  applicable  to  other  quasi-neutral  3heet  configurations 
and  another  example,  f(z)  -  z/5,  i3  discussed  in  C-P. 

We  onoose  the  gauge  such  that  the  vector  potential  is  Ay(x,z)  * 
-33?(z)  ♦  3„x ,  where  dF(z)/dz  •  f(z).  The  3lngle-particle  motion  is 
described  by  the  equation  of  motion 

dv 

-  ^v  x  3  (2) 

dt  c—  — 

This  vector  equation  possesses  three  exact  constants  of  motion;  the 

Hamiltonian  H  -  mv^/2  where  *  v  ^  *  v  ^  *  v  ^ ,  the  canonical  momentum 
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0  and 


However,  for  C.t  and  ?^, 


[C  ,  ?  :  •  -  mC  (3) 
x  y  n 

so  that  Cx  and  ?y  are  not  in  involution.  This  indicates  that  the  particle 
orbits  may  be  stochastic.  As  a  general  remark,  a  Hamiltonian  system  with  N 
degrees  of  freedom  is  ir.tegrable  if  and  only  if  there  exist  N  constants  of 
motion  that  are  in  involution.  Physically,  tr.e  existence  of  such  global 
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cor.3tar.ts  of  motion  means  that  one  can  find  a  canonical  transformation  to  a 
frame  ir.  which  the  M  coordinates  are  cyclic.  The  present  system  is  an 
interesting  example  of  a  nonintegrable  system  which  possesses  three  exact 
constants.  For  more  detailed  properties  of  equation  (2),  see  C-P. 

In  describing  the  orbits,  the  following  normalization  will  be  used;  b«, 

«  3„/3  ,  X  -  (x  -  ?  /mfi  )/b  <5,  Y  -  (y  ♦  C  /tfj/b  J,  Z  -  z/S  5,  t  *  li  t, 

■**  2  ?  2  ^  *  **  " 
and  H  -  H/(mb  ).  A  useful  technique  for  disolayir.g  the  long-time 
n  n  . 

properties  of  the  particle  motion  is  to  use  the  Poincare  surface  of  section 
method  (3ee,  for  example,  Lichter.berg  and  liebermar.,  1983).  For  our 

A 

purpose,  a  surface  of  section  plot  at  2  »  0  for  a  given  value  of  H  is 

constructed  by  following  an  orbit  by  numerical  integration  and  recording 

the  coordinates  X  and  X  *  dX/dn  at  each  point  where  the  orbit  crosses  the 

equatorial  plane.  Figure  1(a),  reproduced  from  C-P,  shows  such  a  plot 

for  H  -  500  and  b„  -  0.1,  evaluated  at  2  •  0.  All  kinematically  allowaole 

*  1/2 

orbit3  are  confined  within  the  circle  of  radius  (2H) 

The  orbit3  in  the  region  marked  A  are  bounded  and  ir.tegrable.  There 
exists  an  additional  invariant  in  this  localized  region  of  the  phase 
space.  This  means  that  the  motion  of  an  i.ntegrable  orbit  is  constrained  to 
a  two-dimensional  invariant  surface  whose  cross-section  through  Z  ■  0  i3  a 
closed  curve.  However,  this  additional  invariant  is  not  a  global  isolating 
constant  and  is  not  expressible  in  closed  form  ir.  terms  of  elementary 

functions.  Mote  also  that  ir.tegrable  orbits  are  not  necessary  adiabatic. 

The  figure  also  clearly  demonstrates  that  there  is  a  large  stochastic 
region,  marked  3.  The  stochastic  region  is  disjoint  from  the  ir.tegrable 
region  A;  there  is  r.o  orbit  that  can  connect  the  two  regions.  The  oroits 
ir.  region  3  are  stochastic  in  the  sense  that  orcit3  are  sensitively 

dependent  on  initial  conditions,  with  two  nearby  orbits  diverging  rapidly 
with  time.  Far  from  tne  equatorial  plane  (z  >>  $),  the  motion  is 
-egular.  We  note  that  Fig.  3  of  Wagner  et  al .  :*?79)  ccr-eapo.nds  it  a 

stochastic  orbit. 


The  nonintegrable  orbits  may  be  thought  of  as  forming  two  flux  tubes, 
which  are  mirror  images  of  eacn  other,  originating  from  arc  escaping  to 


Infinity. 


'he  ~egions  C  inrough 


are  the  1 


0  cross-sections  of  the 


’lux  tubes  as  they  thread  tne  equatorial  plane .  A.l  crcits  at  infinity 


that  car.  reacn  Z  •  0  are  mapped  into  region 


The  orbits  then 


succe33.vely  cros3  regions  Z2  through  15.  These  regions,  11  through  15, 
have  interesting  substructures  as  shown  in  Fig.  i(b).  A  subset  of  the 
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orbits  entering  Cl ,  those  crossing  St  and  71 ,  enter  the  stochastic  region  3 
after  crossing  S5  or  75.  7hat  is,  region  3  can  only  be  accessed  from 
region  C5.  7he  remainder  of  orbits  escape  to  infinity  after  crossing  C5  or 
C4  just  above  74.  7he  latter  orbits  are  referred  to  as  transient  orbits  by 
C-P.  7hese  transient  orbits  appear  to  be  the  type  studied  by  Speiser 
(1967).  Only  orbits  of  this  type  are  shown  in  Figure  1(b).  In  Fig.  2,  we 

A 

3how  a  surface  of  section  plot  for  H  -  500  and  bn»  0.15.  Note  that  all  the 
basic  features  are  similar  to  those  in  Fig.  1(a)  (bn  -  0.1).  However,  with 
a  stronger  3r  component,  a  larger  fraction  of  the  phase  (i.e.,  the  regions 
corresponding  to  Cl  through  C5)  i3  occupied  by  the  flux  tube  structure  and 
the  transient  orbits  therein.  7his  is  because  the  field  lines  are 
"straighter"  for  larger  bn  so  that  the  orbit3  can  translate  in  the  z 
direction  more  easily.  (In  the  limit  b^  •*  •,  all  orbits  move  in  the  z 
direction  freely. )  Similar  structures  exist  in  the  surfaces  of  section  for 

A  A 

different  values  of  H.  In  general,  as  the  value  of  H  i3  reduced,  large- 
scale  integrabie  regions  become  more  fragmented  and  complicated.  For  bn  » 
0.1,  large  integrabie  regions  vanish  for  H  <  6.2  and  essentially  all  orbits 
are  of  the  stochastic  or  the  transient  types.  For  the  parabolic  case, 
similar  features  exist.  A  more  detailed  description  of  the  phase  space 
structures  for  both  cases  can  be  found  in  C-P. 

III.  Current  Distributions 

The  surface  of  section  plots  provide  another  important  piece  of 

information.  It  is  clear  that  the  integrabie  (regular)  orbits  in  region  A 
carry  no  net  current.  On  the  other  hand,  the  nonintegrable  orbits,  i.e., 
the  transient  and  stochastic  orbit3,  do  carry  currents  because  they  come 

A 

from  ar.d  escape  to  infinity.  For  example,  for  H  -  500  in  the  Harris  oase, 
all  nonintegrable  orbits  (ions)  enter  the  equatorial  plane  through  region 
Cl  (Fig.  1(a))  and  escape  from  Co  and  the  nearby  hashed  regions.  This 
means  that  ''egardies3  of  the  actual  path  for  different  orbits,  there  i3  a 
net  drift  in  the  +'t  direction  in  the  vicinity  of  the  equatorial 
plane  ( j Z  j  i  o).  The  electron  motion  is,  of  course,  the  opposite  so  that 

the  electron  and  ion  currents  are  additive.  For  the  parabolic  case  in 

whion  essentially  all  nonintegrable  orbits  are  eventually  reflected  back, 

the  orbits  still  3et  up  a  r.et  drift  near  the  equatorial  plane.  In  this 
case,  the  return  drift  is  established  away  from  tne  C  «  0  plane  with  zero 
total  drift  integrated  over  ail  Z.  Thus,  for  both  the  Harris  ar.d  parabolic 
case,  the  non- integrabie  orbits  form  a  net  current  in  the  dawr.-to-dusk 
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direction  near  the  equatorial  plane.  Note  that  Stern  ar.d  Palmadesao  (1975) 
showed  that  particles  which  are  trapped  between  mirror  points  exhibit  no 
net  drift  across  the  equatorial  plane.  Subsequently,  Cowley  (1973)  and 
Pellat  and  Schmidt  (1979)  generalized  the  proof.  From  our  analysis,  we 
conclude  that  the  proof  is  exactly  applicable  to  the  integrable  orbits. 
The  transient  and  stochastic  orbits  in  the  Harris  case  do  carry  a  net 
drift  because  the  motion  is  unbounded  in  the  z  direction.  In  the  parabolic 
field,  the  proof  is  still  applicable  but  there  is  a  net  drift  near  the  Z  * 
•3  plane  as  described  above. 

17.  Generation  of  Free-Er.ergy 

A  prominent  feature  of  the  particle  motion  in  the  quasi-neutral  sheet 
geometry  is  the  existence  of  disjoint  regions  in  the  phase  space  near  Z  ■ 
3,  each  region  consisting  of  orbits  of  distinct  nature;  (1)  bounded 
integrable  orbits,  (2)  'unbounded  stochastic  orbits  and  (3)  unbounded 
transient  orbits.  In  effect,  there  exi3t  boundary  surfaces  between  the 
disjoint  regions  which  do  not  allow  randomization  of  information  or 
particle  energy  because  of  the  .phenomenon  of  differential  memory  (Chen  and 
Palmadesao,  1985).  If  noise  fields  are  introduced  to  the  magnetic  field, 
the  boundary  surfaces  will  break  up  and  allow  ’’diffusion"  of  orbits.  The 
time  scale  for  this  process  depends  on  the  strength  of  the  noise  field  ar.d 
13  slow  for  low  levels  of  noise.  The  basic  disjointness  of  the  regions  and 
differential  memory  are  expected  to  be  insensitive  to  low-level  noises. 

The  concept  of  differential  memory  has  a  number  of  physically 
interesting  implications.  Suppose  the  system  contains  a  population  of 
charged  particles  in  thermal  equilibrium.  There  must  also  be  a  steady 
supply  of  particles  in  the  distant  regions  to  maintain  equilibrium.  If  the 
oaraaeters  of  the  distant  plasma  distribution  are  changed,  then  the  fact 
the  oroit3  are  divided  into  di3tl.nct  types  belonging  to  disjoint  -esions  of 
the  pr.ase  3pace  -esuits  in  a  hignly  non-Maxwell  tan  distribution  oecause  of 
differential  memory.  The  resulting  ncn-Maxwellian  distr ibutions  may  oe 
written  in  the  fora  F(H,?y)  for  the  normtegrafcle  regions  ar.d  FCH.P^.g)  for 
the  integrable  regions  where  g  is  the  additional  invariant  of  motion ,  The 
existence  of  boundary  surfaces  car.  thus  establish  free-er.ergy  wnich  car  be 
tapped  by  suitable  plasma  instabilities. 

In  the  context  of  the  earth’3  magr.etotail ,  it  is  believed  that  t re 
coil isionless  tearing  instability  may  play  an  important  role  in 


reconnection  processes.  It  has  recently  been  shown  (Chen  and  Palaadesso, 
1984b;  Chen  and  Lee,  1985)  that  the  growth  rate  of  the  collisionless 
tearing  node  in  a  non-Maxwell i an  neutral  sheet  can  be  substantially 
greater,  by  a  few  orders  of  magnitude,  than  in  the  Maxwellian  case.  The 
general  fora  of  the  distribution  functions  which  can  provide  the  necessary 
free-er.ergy  is  ?(H,Py,C)  where  C  is  an  independent  constant  of  notion. 
Thus,  if  the  aagr.etotail  which  is  initially  in  theraal  equilibrium  is 
subjected  to  changes  in  external  conditions,  e.g.,  the  3olar  wind,  pressure 
distribution,  zagr.etopause,  etc.,  then  the  magnetic  topology  can  develop 
appropriate  plasma  distributions  to  drive  the  non-Maxwellian  tearing 
mode.  An  important  point  is  that  the  free-energy  for  this  instability  is 
the  particle  energy  and  that  the  conventional  tearing  mode  also  occurs  but 
on  faster  time  scales.  That  is,  changes  in  external  conditions  can 
generate  by  the  process  of  differential  memory  an  additional  source  of 
free-energy  in  the  form  of  non-Maxwellian  distributions. 

We  add,  however,  that  the  above  discussion  does  not  take  into  account 
the  influences  of  the  normal  magnetic  field  on  the  instability  itself  under 
certain  -circumstances  (Galeev  and  Zelenyi,  1976;  Lembege  and  Pellat,  1982; 
Coroniti,  I960).  This  point  warrants  further  consideration  and  will  be 
addressed  in  a  separate  paper.  Finally,  the  above  results  and  discussions 
are  based  on  3 ingle- particle  motion.  In  order  to  determine  more 
quantitatively  the  nature  and  associated  time  scales  of  these  processes,  it 
is  necessary  to  follow  an  ensemble  of  particles.  Thi3  work  is  currently 
underway. 
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